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1 INTRODUCTION 



Since being first predicted by iBaade fc Zwickvl l|l934 ). 
neutron stars have been observed to display a wealth 
of phenomena. Over the past decades a new class 
of neutron stars has emerged, through the studies of 
the emission mechanism of soft 7-ray repeaters (SGRs) 
and anomalous X-ray pulsars (AXPs). Termed mag- 
netars, they are powered by their extreme magnetic 
fields, typically of the order 10^*-10^^ G, rather than 
their spin-down energy loss as in the case of pulsars 
(■Paczvnski' '1992'; 'Usov '1992'; 'Duncan fc Thompson! 1 19921 : 
[Thompson & Duncan 1993, 1995, 1996) (for an overview , 
see, e.g.. [Thompson fc Duncanll200ll : iHarding fc Laill2006l ). 
Now strengthened by observational measurements in- 
cluding their slow spin periods, fast spin-down rates 
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1998"; 'Woods & Thompson' '2006') and 
spectral propert ies I Gavriil ot al. 2002; Ibrahim ct al. 20031: 
iRea et al.ll2003l ). there is increasing evidence for their ex- 
treme magnetic fields (for a review of observations, see, e.g., 
I Woods fc Thompson|[200^ . 

The origin of the magnetic field however remains de- 
bated. Of the magnetar candidates — 5 SGRs, 6 AXPs, and 
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ABSTRACT 

Under the fossil field hypothesis of the origin of magnetar magnetic fields, the magnetar 
inherits its magnetic field from its progenitor. We show that during the supernova of 
such a progenitor, protons may be accelerated to '--^lO* GeV as the supernova shock 
propagates in the magnetic stellar envelope. Inelastic nuclear collisions of these protons 
produce a flash of high-energy neutrinos arriving a few hours after thermal (10 MeV) 
neutrinos. The neutrino flash is characterized by energies up to O(IOO) GeV and 
durations seconds to hours, depending on the progenitor: those from smaller Type Ibc 
progenitors are typically shorter in duration and reach higher energies compared to 
those from larger Type II progenitors. A Galactic Type lb supernova leaving behind 
a magnetar remnant will yield up to ~160 neutrino induced muon events in Super- 
Kamiokande, and up to '^7000 in a km'^ class detector such as IceCube, providing a 
means of probing supernova models and the presence of strong magnetic fields in the 
stellar envelope. 

Key words: acceleration of particles - neutrinos - stars:magnetic fields - super- 
novae:general - pulsars: general 



a few radio pulsars — approximately a third are associated 
within known young supernova remnants, suggestive of an 
origin in ma ssive star explosions ( Woods fc Thompsonll2006r ) 
(but see also lGaensler et al.ll2001 ). The ages of the remnants 
are ^10* years, consistent with the inferred ages of magne- 
tars derived from their spin-down rates ("spin-down age"), 
which lie tightly between 10^ and 10^ years. If magnetars 
are young neutron stars resulting from core-collapse super- 
novae, then their magnetic fields could have been inherited 
from their progenitors. This is the so-called fossil field hy- 
pothesis, where magnetic flux is conserved and the field is 
ampl i fied during the core-collapse proc ess (JHarding fc Lail 
120061 : iFerrario fc Wickramasinghd I2OO61 ) . Observations re- 
veal that there is, in principle, enough magnetic fiux, i.e., 
the magnetic fiux of O- stars, derived from their recen tly de- 
tecte d magnetic fields JDonati et al.ll2002l . [2OO6 ; Petit et all 
l200a ). are comparable to those of magnetars. On the other 
hand, the magnetic field may be generated by a convective 
dyna mo in the first 0(10) seconds of the protoneutron star 
birth (IThompson fc Duncanlll9 93r). In this process, the en- 
ergy in differential rotation is converted to magnetic energy. 
It is not yet clear which is the dominant process. 

Here we discuss proton acceleration and production of 
high-energy neutrinos during a supernova which leaves be- 
hind a magnetar. We focus on the fossil field hypothesis, 
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which endows the star with strong stellar magnetic fields. 
We explore three progenitors corresponding to Type Ic, lb 
and II supernovae, and show that proton acceleration is re- 
alised in all the cases during the propagation of the shock 
through the stellar envelope. The maximum proton energy 
is sufficiently high that neutrinos are produced through in- 
elastic proton-proton collisions. The neutrinos have energies 
of 0.1 GeV up to O(IOO) GeV, which is significantly higher 
than thermal neutrinos from the collapsed core (~10 MeV). 
The signal thus resides in a unique energy window, a positive 
point from a detection perspective. 

The neutrino emission and its detectability depends on 
the energy loss rate of muons and pions. We find that a su- 
pernova explosion at the Galactic Centre results in 70-160 
neutrino induced muon events in Super-Kamiokande, with 
the largest value for a Type lb supernova. Higher energy neu- 
trinos in the sensitivity range of km'' class detectors are more 
strongly dependent on the progenitor, with 200-6600 events 
per supernova. The progenitor dependency is the product of 
two physical processes, proton acceleration and energy loss 
of pions and muons. For Type Ic supernovae, the compact 
progenitor results in strong cooling of pions and muons, so 
that high-energy neutrino emission is suppressed. The cool- 
ing is dominated by inverse-Compton scattering on electron 
synchrotron photons. For Type II supernovae, the maximum 
neutrino energy is intrinsically low. The greatest emission of 
high-energy neutrinos is realised for a Type lb supernova. 

In all the cases, the neutrinos arrive tens of hours af- 
ter thermal neutrinos, and last between seconds to hours 
depending on the progenitor radius. The background atmo- 
spheric neutrino is significantly smaller, and detection is es- 
sentially background free in most cases. Given the nature 
of the fossil field hypothesis, detection of these high-energy 
neutrinos can provide support for the stellar origin of mag- 
netar magnetic fields. Detection also provides useful diag- 
nostics for supernova properties. 

The paper is structured as follows. In section 2 we dis- 
cuss the progenitor density profile and the supernova shock 
environment. These are the required background for inves- 
tigating proton acceleration in the shock, which we address 
in section 3. We also include the production of neutrinos 
and their detection prospects in section 3. In section 4 we 
finish with discussions. Throughout this paper we define 
Qa ~ Q/10" for a quantity Q in cgs units. The exception is 
Ej;, the energy of particle species x, for which we use GeV. 



2 SETUP 

In this paper we investigate proton acceleration at a super- 
nova shock propagating through the stellar envelope. Both 
proton acceleration and particle cooling depend on quanti- 
ties near the shock front, such as the magnetic field, pho- 
ton density, and particle density. We therefore start with 
a discussion of these required quantities, and treat proton 
acceleration in the next section. Here, we discuss the pro- 
genitor density profile (section 2.1), the progenitor magnetic 
field (section 2.2), and the particle density and temperature 
around the supernova shock (section 2.3). 
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Figure 1. Plot showing "acceleration possible" (below dotted 
lines) and "no acceleration" (above dotted lines) regions, where 
the numbers by the dotted lines indicate the maximum acceler- 
ated proton energy. The density profiles of our chosen models 
are also plotted: CO star (resulting in Type Ic supernova, dot- 
dashed). He star (resulting in Type lb supernova, dashed), and 
BSG (resulting in Type II supernova, solid). The supernova shock 
traces the density curves. Therefore, acceleration becomes possi- 
ble as the shock approaches the stellar surface. 



2.1 Stellar density 

First we discuss our selection of progenitor density profiles. 
It is widely believed that neutron stars are formed from 
the collapse of massive OB stars, with main-sequence mass 
ranging between 8 < M/Mq < 45. We assume that mag- 
netars are similarly produced as remnants of core-collapse 
supernova of massive stars. Due to the lack of observational 
evidence for which type of core-collapse preferentially pro- 
duces a magnetar remnant [j we consider three supernova 
progenitors: a CO star, a He star, and a blue supergiant 
(BSG), which respectively give rise to a Type Ic, lb and II 
supernova. 

For the CO star model (Type I c progenitor) , we 
adopt the presupernova model 16SK from lWooslev fc Hegerl 
l|2006t ). motivated by the 7-ray burst association with core- 
collapse supernovae. This is a rapidly rotating star with solar 
metallicity and initial mass of 16 Mq. Due to strong mass 
loss, the final mass is Af, ~ 5 Mq and radius is i?, « 5 x 10^" 
cm. The outer envelope is radiative and dominated by car- 
bon. 

For the He star model (Type lb progenitor) , we 
adopt the presupernova model 12SE from lWooslev fc Hegerl 
(200u). This is a rapidly rotating star with solar metallicity 
and initial mass of 12 Mq . Due to mass loss, the final mass 



^ However, we note that a connection between magnetar birth 
and type Ic supernova h as been proposed in the context of lo w- 
luminosity 7-ray bursts l lMazzaU et al.ll200^ : iToma et al.ll2007l) . 
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is M, ^ 7Mq and radius is i?* Ri 2 x lO" cm. The outer 
envelope is radiative and dominated by helium. 

For the BSG model (Type II progenitor) we 
adopt the non-rotating presupernova model 16TA from 
IWooslev fc Hegeij (120061 ). This is a low metallicity star with 
an initial mass of 16 Mq . There is little mass loss, and the 
final mass is M* « 16 Mq and radius is i?, ~ 3 x lO^'^ cm. 
The outer envelope is radiative and dominated by hydrogen. 
Note that we only consider BSGs which have radiative en- 
velopes, and do not consider more common red supergiants. 
This is because radiative envelopes are necessary for su rvival 
of fossil magnetic fields (|Mossll2003l : iTout et al.ll2004l ). 

In Fig. [l] we show the density profiles of our adopted 
models. Throughout this paper, when we provide numerical 
values to equations, we adopt the He star model and the ra- 
dius 0.9i?, , for illustrative purposes. The illustrative density 
is ~10^''' gcm~^. 



2.2 Stellar magnetic field 

Now we discuss the progenitor magnetic field. We work un- 
der the fossil origin for the magnetic field of neutron stars. As 
we discuss below, this requires the progenitor to be strongly 
magnetized. 

The magnetic fields of massive stars have recently been 
directly detected. The measured strengths of surface dipole 
fields are as high as 1 kG. Although at the present time there 
are only a handful of such direct detection, including 9 Ori 
C with 1.1 ±0.1 kG JPori ati et al. 2002), HD191612 with 
1.5 kG (|Donati et al.l 120061 ) . and two in the Orion Nebula 
with llOOl^oo G and 650ti7g G (iPetit et al.ll2008h . it has 
been speculat ed that magnetism may be widespread among 
massive stars (jPetit et al.ll2008r ). Observed properties of the 
magnetic field, such as its global nature and no clear corre- 
lation with stellar parameters, favour a fossil origin over a 
dynamo origin; i.e., the magnetic fields are fossil remnants 
from the star formation stage, as relics of the field that per- 
vaded the interstellar medium (|Donati et al.ll2006l '). 

Interestingly, the magnetic flux of 9^ Ori C, calcu- 
lated from its observed magnetic field and inferred radius 
ISimon-Diaz et al.ll200^ to be (7 ± 3) x lO^'^Gcm^ is re- 
markably close to that of the highest field magnetar SGR 
1806-20, ~3 X 10^* G cm^ (assuming a radius 10*^ cm) . There- 
fore, in principle there is enough magnetic flux present in a 
massive magnetized star to explain the magnetic fields of 
magnetars, and hints at a possible evolutionary link. This is 
the so-called fossil field hypothesis of neutron stars, where 
magnetic flux conservation results in field amplification dur- 
ing collapse. In addition to the assumption that magnetic 
flux is conserved, it is implied that the fleld must some- 
how survive the post-main-sequence evolution and the var- 
ious internal structural changes during the formation of 
the neutron star. Despite these ideal assumptions, the sce- 
nario could provide a powerful explanation fo r the wide 
range of magnetic flelds present in neutron stars (IMosJuOOSl : 
iTout et al.ll2004l : lFerrario fc Wickramasinghell2005l . l2006l ). 

It is worth commenting that the fossil fleld hypothe- 
sis of neutron stars is an extension of the fossil theory for 
magnetic white dwarfs. In the case of white dwarfs, the frac- 
tion of the strongly magnetized population is roughly com- 
patible with the fraction of m agnetism in early-type stars 
from which white dwarfs arise l|Wickramasinghe fc Ferrarid 



|2005|). The broad study of magnetic flelds in massive OB 
stars has recently started (jPetit et al.l 120081 ) . which would 
enable similar studies for neutron stars. 

Under the fossil field hypothesis, the neutron star inher- 
its its magnetic field from the progenitor. First we define a 
mass cut, Mcut < 2M0, corresponding to the requirement 
that the collapsed object forms a neutron star and not a 
black hole, i.e., Mcut is the collapsed mass, while the overlay- 
ing mass M* — Mcut is ejected by the supernova shock. The 
magnetic field of the remnant object must therefore originate 
from the material within the mass Mcut • In our adopted pro- 
genitor models, the iron core mass is Mpe « 1.5 Mq, which 
approximately corresponds to the required mass cut. We 
therefore assume A/cut ~ A/pc for simplicity (the exact lo- 
cation is not critical for our purposes). Thus, the magnetic 
field at the surface of the iron core is found to be, from 
conservation of magnetic flux, 

where _Rfo ~ a few xlO* cm is the radius of the iron core. 

The magnetic field in the stellar interior, including the 
envelope, is largely unknown. Following this difficulty, we pa- 
rameterize the field strength according to a power-law with 
index n, 



B{r) = Be 



Rfc 



G for _Rfo < r < R,, 



(2) 



where we have taken the normalization at the iron core ra- 
dius. For a dipole field, n = 3. In the fossil theory for neutron 
stars that we adopt throughout the preset paper, n = 2 is 
obtained by equating Eq. ([2]) to the magnetic field strength 
measured at the stellar surface. We note that for our illus- 
trative model, the He star, the resulting magnetic field is 
^W^-^ G at radius 0.97?,. 



2.3 Shock environment 

As the shock propagates through the star, it accumulates 
stellar material, shocking it. Across the shock, material is 
compressed; the downstream (labelled d) density is pd ~ 
Kpu, where pu is the stellar (upstream) density discussed in 
section 2.1, and n is the compression ratio. We adopt the 
fiducial value of k = 4, but note that in principle k can 
reach up to 7 for a radiation shock. 

For compact stars, the supernova shock can reach 
velo cities of ^O.lc as it leaves the stellar surface 
( Mat zner fc McKed 1 19991 ) . Similar velocities are also seen 
in numerical studies of explosions of O-Ne-Mg cores, where 
the steep densi ty gradient drives a f ast shock velocity as high 
as 10^° cms~^ IjKitaura et al.ll2006l ). Since steep density gra- 
dients are generically expected near surfaces of stars, in the 
current paper we adopt Vs = 10^'^ cms~^. 

In the optically thick interior of the p rogenitor, the 
supernova shock is dominated by radiation (|Weaveiill976l : 
iMatzner fc McKeel ll999). The radiation temperature is ap- 
proximately given by the relation aT^ ~ puV^/2, where 
a = 4ct/c and a is the Stefan-Boltzmann constant, so that 



Tr ~ 800p;' 



1/4 1/2 



eV. 



(3) 



The energy densities of the photon and magnetic fields are 
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given as Ur ~ pu'ul/'i and Ub ~ B^ /Sir, respectively. Quite 
generally, Ur S> Ub', this implies the relative importance of 
inverse-Compton over synchrotron process. 



3 NEUTRINOS FROM SUPERNOVA SHOCKS 

3.1 Proton acceleration 

In this section we discuss proton acceleration by diffusive 
shock acceleration, which produces a power-law spectrum 
dN/dp oc p"" for relativ i stic protons, where p is the par- 
ticle momentum dOrurvl 11983 : iBlandford fc Eichled 1 19871 : 
iMalkov fc Drurvl I2OOII ). The acceleration time is tacc ~ 
lOD/v'^, where D = r/rg/S is the diffusion coefficient of par- 



ry is a factor which depends on the ratio of energy in the 
ordered magnetic field to that in the turbulent magnetic 
field. Although there is some uncertainty regarding this co- 
efficient, due to the high level of turbulence close to the 
shock front, we assume diffusion in the Bohm limit (r; — > 1). 
We note that this is the most efficient limit for acceleration. 
We then obtain (in the relativistic regime) 

10 ro 



fa 



l3_ 

3 i;? 



1X10 ep,iW^^g.5B4.5S, 



(4) 



for our fiducial parameters, where ep is the proton energy 
in GeV (we reserve Ep for the total energy in accelerated 
protons). In the plasma, the proton acquires energy from 
the bulk plasma motion through resonances (jLongair 1994). 
The plasma waves taking part in this interaction are caused 
by oscillating electrons. If the photon-electron collision rate 
exceeds the plasma frequency, these waves will not be pro- 
duced. Fortunately, we can show that in the radiation domi- 
nated shock near the stellar surface, the photon-electron col- 
lision rate is smaller than the plasma frequency: while the 



photon-electron collision rate is Vj 



CUjO-T = 3 X 10* s ^, 



the plasma frequency is Upi — 2 x 10 s~ . Here, ar is the 
Thomson scattering cross section and n-y ~ 0.33(Ur/ch)'^''^ 
is the photon density. 



1 ,, iri22 3/4 3/2 

1 X 10 P,'rV'Q. 



(5) 



Proton acceleration must compete against the shortest 
of three time scales, namely i) proton escape from the ac- 
celerating region, ii) proton energy loss time, and Hi) age of 
the accelerator. In supernova shocks, acceleration is limited 
by the proton energy loss time. 

Accelerated protons lose energy by interacting with the 
dense photon field and dense protons in the shock vicin- 
ity. We first discuss proton cooling by inelastic collisions 
with target protons. We adopt the target proton density 
Kpu/nip {nip is the proton mass), the prot on-proton coUi- 
sional cross section app « 4 x 10~^^ cm^ (|Eidelman et al.l 
[2OOJ), a nd that 20% of the proton energy is lost in each 
collision (|Bellandi et al.lll994h . The cooling time, defined as 
ep/|dep/df|, is then 



^»» — 



0.2cappKpu/mp 



2xlO-X^_5S. 



(6) 



Protons may also lose energy by pair-production on the 
photon field, also known as the Bethe-Heitler (BH) process. 
The pairs are produced at rest in the centre of mass frame, 
and the energy lost by the proton is Atp = 2meC?^cm, 



where 7cm = (ep + e^)/{rnlc^ + 2epe^)^''^. The energy 
loss rate is then dep/df 



f[(28/9)ln(2epe.y/mpr 



n^cuBH^f-p, where ubh = 
— 106/9] is the cross section. 



Here, a and r^ are the fine structure constant and classical 
electron radius. BH cooling occurs only for protons above 



5 X io>; 



-1/4 -1/2 



GeV. Its time scale, determined 



assuming the average photon energy e^ ~ 2.7Tr and substi- 
tuting Eq. (O, reaches a minimum near ep ~ 10^ GeV, but 
is slower than the coUisional cooling time. 

Finally, protons may cool by synchrotron and inverse- 
Compton processes with the magnetic and photon fields, 
respectively. However they are slow, with cooling times of 



3m„c' 



4„3 



fic = 



4(TTm|ep[/s 
Sm-pC^ 
AaTTn'iepUr 



4 X lO^e'lB. 



2 
4.5 s. 



'400e„.iP„_^_5U.g.5S, 



(7) 



(8) 



respectively {m^ is the electron mass). 

We conclude from the above analysis that proton- 
proton coUisional cooling dominates the cooling of protons. 
Equating this cooling time to the acceleration time yields 
the maximum proton energy. 



vhB 



0.2cappKpu/r. 



2 4 
rripC 



1/2 



(9) 



which is ep^max ~ 160 GeV for our fiducial He star and ra- 
dius 0.97?,. Note that once the shock leaves the star, the 
proton density is expected to fall dramatically, and cooling 
will most likely be dominated by pair-production and radia- 
tion processes. Inside the star however, Eq. ^ may be safely 
used. From equation Eq. ^ one can identify "acceleration 
possible" regions and "no acceleration" regions on the B-pu 
plane. We show this in Fig. [ij making use of Eq. ([2]) to plot 
the radius on the horizontal axis. We plot three threshold 



lines corresponding to tp 



= 1.23, 10, and 100 GeV. We 



also show the stellar density profiles adopted. We see that 
proton acceleration becomes possible as the shock nears the 
stellar surface. 

Finally, we add a note regarding the possible effects 
of shock-accelerated electron synchrotron photons. Given 
the large magnetic field, the energy in accelerated electrons 
will be rapidly converted to synchrotron photons. An in- 
crease in the photon density can potentially cause two prob- 
lems. First, the photon-electron collision rate may exceed 
the plasma frequency. Second, proton cooling by inverse- 
Compton will become faster. We find that for typical energy 
fractions of relativistic electrons (less than a few percent of 
the shock energy fivp (lAllen et al.ll200ll : iBamba et alll2003l : 
IWaxman fc Loebll2001^ '). the increase in photon density is 
not sufficient to cause these problems. For example, sub- 
stituting Ur = ^eEaxp/V where V is volume, and conser- 
vatively assuming ^e = 0.01, the inverse-Compton cooling 
time scale time is 



flC.e ~ 60ep_ii?,_ii,3^c,-2-£'Gxp,5lS, 



(10) 



which is much slower than coUisional cooling. Moreover, 
these conclusions remains unchanged even if the energy frac- 
tion in electrons is maximised (i.e., the unrealistic case where 
S,e = 1). We therefore safely ignore this effect for protons. 
Note that it is however important for pions and muons, 
which we discuss below. 
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3.2 Neutrino production 

Inelastic proton-proton interactions producing pions occur 
for protons witli energies above tiie tliresiiold energy ep,th = 
{^{mjy + nrin + ?n^)^ — m^/rap = 1.23 GeV. Tiie decay of 
ciiarged pions produce neutrinos tiirougii tt"^ — + /i"*^!^^ and 
7r~ -^ fi~i>fj,. Tlie muon neutrino energy is ei, ~ 0.25e,r ~ 
0.05ep, since e^ ^ 0.2ep. Competing witli this decay is 
pion cooling by radiative and collisional processes. We de- 
fine them similarly as we did for protons , but by adopt- 
ing cr^p = 5 X 10~^®cm^ JEidelman et all |200 j ) and that 
80% of the pion ener gy is lost in each pion-proton collision 
(jBrenner et al.lll982l ). We find that radiative cooling domi- 
nates, in particular inverse-Compton scattering on electron 
synchrotron photons. Equating the dominant cooling time 
to the decay time, r,r ~ 4 x 10~^ep,i s, we define the break 
energy e\^^l- Below this energy, pions decay without energy 
loss, while above this energy the decay spectrum is sup- 



pressed by a factor ijr.rad/Tir 



To take this spectral 



steepening into account, we define the suppression factor 



C(^.) 



1*= I/, brk/ *'''-' 



for e^ < e 
for e^ ^ e 



iy.brk 

(,r) 



(11) 



where e 



(vr) 



300 GeV for our chosen parameters and He 
star (f»40 GeV for the CO star and f«2 x lO" GeV for the 
BSG). 

In high-energy proton-proton interactions, tt^, tt", and 
tt" are produced in nearly equal numbers. While inside the 
star, the density is high enough for the proton-proton optical 
depth Tpp = /fl * drapppu{r)/mp to be very high, so that all 
of the energy in accelerated protons is converted to pions. 
Thus, we normalize the pion spectrum by the total energy, 
rather than by the particle number. Since 2/3 of the pro- 
duced pions are charged, the total energy in charged pions 
is E^± « 2Ep/3. Here, Ep is the total energy in accelerated 
protons, which we parameterize as Ep = ^p-Eoxp, where ^p is 
the fraction of the supernova shock energy -Eexp that is chan- 
nelled into accelerated protons. For strong shocks around su- 
pernov a remnants, the inferred va l ues of ^p are of the order 
of 0.1 (|Blandford fc Eichlerl HqStI : IffiUad ^OOSl. In charged 
pion decay, the muon neutrino takes approximately 1/4 of 
the pion energy. The flavour ratio of neutrinos produced at 
the source is u'^ : u^ : ly'^ — : 1 : 0. Neutrino oscillations 
en route to a detector on Earth then lead to an observed 
ratio Ve '■ I'n '■ i^T ~ 1 : 1-8 : 1.8 (jLearned fc Pa kvasa 1995; 
iBeacom et al.ll2003l : ?). The differential fluence of muon neu- 
trinos {i'i_i + i'fi) from a supernova at a distance D is therefore 
given as 

de, ^ 4nD^ln{ep,rn../ep,ti.)el^^'"'' ^ ' 

where C(^i') is the suppression factor due to pion cooling, 
^^ = 1/15, and ^th is the fraction of Ep lying above the pion 
production threshold energy ep,th, 



Cth — 



j;-'""" depepdN/dp 



r 

Jrr. 



dep€pdN/dp 



(13) 



Muon decay also contributes to the neutrino flux, 
through /!+ ^ e^VfiUe and /i^ -^ e^v^Ue- However, due 
to their smaller mass and longer decay time, they only sig- 
nificantly contribute at lower energies compared to pions. 



The break energy, determined in the same way as for pi- 
ons, is fj'^l-^y^ ~ 20 GeV for our chosen parameters and 
He star («3 GeV for the CO star and «1300 GeV for the 
BSG). Note that m uon collisional cooling, eva luated using 
the cross section of iBorog fc PetrukhinI (|l975l ), is not sig- 
nificant, and muons cool most rapidly by inverse-Compton 
on electron synchrotron photons. Below the break energy, 
the total flavour ratio of neutrinos at the source, i.e., com- 
bined with muon neutrinos from pion decay, is z/° : z/° : 
i/° = 1 : 2 : 0. The observe d ratio after osc illations is then 
I'e : i^/j : i^T ~ 1 : 1 : 1 (|Learned fc Pakv asa 1995), and 
£,v = 1/6. The change in flavour ratio ca n also be a probe of 
this transition (IKashti fc Waxmanll2005h . For our purposes, 
the inclusion of muon decay increases the muon neutrino 
flux by a factor 5/2, for energies below f]^l^^]^- 



3.3 Neutrino-induced muon detection 

In this section we discuss the detection prospects of the neu- 
trinos discussed in previous sections. We address the ques- 
tion of the maximum accelerated proton energy, which deter- 
mines the maximum neutrino energy and hence detectability 
at Super-Kamiokande and IceCube detectors. In particular, 
detec tion by IceCub e is great improved by the IceCube deep- 
core (|Resconill2008l ). with its low detection threshold energy 
of e^ ~10 GeV. 

As the supernova shock propagates outwards, the max- 
imum proton energy increases. In order to address ac- 
celeration in the outer envelo pe, we make use of the 
analytic formula according to iMatzner fc McKed (119991 ) 
for t he mass density at the e dge of the star (see also 
e.g., IWaxman fc Meszarod 120031 '). 



P = P* 



R, 



(14) 



where n = (7 — 1)~^=3 for a radiative envelope, and 7 
is the adiabatic index. We fit this function to our chosen 
stellar models. The CO star is approximately described by 
the parameters p* = 2 gem"'' and R^ — 5 x 10^^° cm. The 
He star is described by the parameters p* = 0.01 g cm~^ 
and Rt = 2 X W^^ cm, while the BSG is described by the 
parameters p, = 6 x 10~^ gcm~^ and 7?* = 3.4 x 10^^ cm. 

Using the fits to the stellar density, we determine the 
maximum proton energy at Rt{tpp), the radius defined by 
a proton- proton optical depth to the stellar surface of Tpp. 
We use Eq. (|9]), which is justified since proton-proton colli- 
sion is the dominant proton energy loss process. This yields 
Ep.max ~ 8 X 10'^ GeV for the fiducial He star at 7?t(t"pp = 5) 
(~1 X 10** GeV for the CO star and ~900 GeV for the BSG). 
Although larger radii yield larger ep,max, we do not consider 
this for several reasons. First, the fraction of accelerated pro- 
tons interacting with target protons rapidly decrease since 
the target density falls as a steep function of radius. Al- 
though protons can also interact with photons, the required 
photon energy is very high, ep£j > 0.3 GeV^. Second, we 
have normalized the neutrino spectrum under the condition 
that all accelerated protons lose energy by multiple proton- 
proton collisions. Third, higher energy neutrino emission is 
typically strongly suppressed due to pion cooling. Therefore, 
while larger radii yield higher energy neutrinos, they do not 
lead to more detected neutrinos. 

The total number of I'^j-induced muon events in 
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a neutrino detector is tiie integral over energy of 
ptVdct<7ccdF,^/de^, where pt is the target density, Vdet is the 
detector volume, and ace is the inelastic neutrino-nucleon 
cross section. We t ake into account the muon range (e.g., 
iBeacom et al.l 120031 ). which effectively increases the detec- 
tor volume by detecting muons produced outside the in- 
strumented volum e. We adopt the cross section used in 
lAshie et al.l (|2005|), which covers the energy range down 
to e^ « 0.1 GeV. Between 1 < e^,GoV < 10^, the cross 
section increases roughly proportional to e^ IjGandhi et al.l 
[1998 ). The energy of the muon faithfully represents the neu- 
trino energy, and we use the average of the f -N an d v^-N 
values of {y) as computed in iGandhi et al.l (|l996l ). where 
y(ey) = 1 — e^/ei/. The opacity of the Earth becomes compa- 
rable to 1 when the neutrino energies are ~10 TeV or more; 
the Earth is therefore totally transparent for the neutrinos 
of our interest. 

We adopt -Ecxp ~ 10^^ erg and distance D = 10 kpc, 
i.e., a supernova occurring near the galactic centre. First 
we discuss detection with Super-Kamiokande. Integrating 
the fluence over neutrino energy O.l-O.OSep^max GeV, the 
expected number of neutrino induced muons is 



Na 



160^p,-l-Ecxp,51-D2 



(15) 



for our fiducial He star. Prospects for detection of the CO 
star and BSG are similar, with total events of 70 and 130 
respectively. Note that the contribution from muon decay 
is non-negligible at Super-Kamiokande energies, since the 
muon break energy is above the detection threshold. For 
example, of the 160 events, we find that ~50 are from muon 
decays. 

For IceCube deepcore, we conservatively estimate the 
effective area as 

2.0 X 10-^(e^/10)^'^ 



Aadie^) = 



6.3(e./102)i 



cm^ 10<e. s^lO^ ,_, 
2 in2 ^ , ^ in4 UOJ 



10^ < e. s; 10^ 



which contains neutrino interaction probability, muo n prop- 
agation, detector response, and event selection (e.g.,|Desiati| 
l2006l : lRottll2008l '). The neutrino energy is in GeV. We inte- 
grate over neutrino energy 10-0. 05ep, max GeV, yielding the 
expected number of neutrino induced muons 



Nu 



66005p,-l£oxp,5ll>2 



(17) 



for the fiducial He star. Total events for the CO star and 
BSG are ~600 and ~200 respectively. The suppression in 
the CO star case is due to the low ej, ^rk' while for the BSG 



it is the low e^ 



Note that we have neglected contribu- 



tions from muon decays for IceCube predictions, which are 
expected to be small. We also note that the break energy is 
a function of ^e, the fraction of total energy in electrons, by 



A^) 



-1/2 



"=i^ brk °^ ^e ■ We have conservatively assumed ^e = 0.01. 
For e xample, for ^e = 10~^ (JAllen et a l. 2001; iBamba et all 
120031 ). the break energy for the Type Ic case is r;10* GeV, 
and the total event number increases from ~200 to ~3000. 
All the events discussed will cluster in a time window 
Ri R,/vs which is of the order of seconds to hours, depend- 
ing on the progenitor size. In comparison, the background 
muon rate due to atmospheric neutrinos, over the entire 
27r steradian sky in 1 day, is ~10 for Super-Kamiokande 
and ~100 for IceCube, using the same assumptions of 
detection efficiency as our signal calculations. These are 
comparable to values in the literature, ~10 for Super- 



Kamio kande (iKaiita fc Totsukal 120011 ) and O(IOO) for Ice- 
Cube IjGonzalez-Garcia et al.ll2005l ). We add that the true 
number of signal events will differ by more than a factor 
of a few, depending on the quality of cuts used. However, 
any increase or decrease from our estimates will similarly in- 
crease or decrease atmospheric background events, and our 
conclusion that the signal overwhelms background is unaf- 
fected. Co mbined with the angula r correlation with an shock 
breakout (jSoderberg et al.l 120081 ) or the optical supernova, 
an essentially background free detection of the neutrino sig- 
nal is possible, as long as A^^ > 1. 

3.4 Diffuse nonthermal neutrino background 

From Eq. (I12p . the emitted number (per unit energy range) 
of nonthermal neutrinos from each supernova is dNv/de^ — 
inD^dF^/de^ ^ 1 x 10^^C(e^,o)Co GeVS adopting fidu- 
cial parameter values for the He star model. These neutri- 
nos have been continuously injected to the Universe since 
the beginning of star formation. Since they are unabsorbed, 
they form a diffuse nonthermal neutrino background in the 
present Universe. Thus, it is of interest to estimate the in- 
tensity of this diffuse background, and compare it with other 
diffuse components such as atmospheric neutrinos and dif- 
fuse thermal neutrino background f rom supernovae (e.g., 
lAndo et al.ll2003l : lAndo fc Satdl2004l ). 

A key element for such an estimate is the occurrence 
rate of supernovae leaving behind strongly magnetized com- 
pact objects in the Galaxy. The relatively young age of 
the observed magnetars indicate a lower Galactic magne- 
tar birth rate of 1 in 10^ years IjKouveliotou et al.l 11994 : 
Ivan Paradiis et al.lll995l ). However, due to low detection ef- 
ficiencies caused by e.g., on/off states of magnetars, it has 
been suggested that the rate may be an order larger, ap- 
proac hing the rate of radio pulsars (jWoods fc ThompsonI 
120061 ). Therefore, we adopt the magnetar birth rate of 
10~''yr~'^ per galaxy. Given the a verage galaxy num- 
ber density, rtg ai ^ 10"^ Mpc"^ (JBlanton et al.l 1200 ll : 

iNakamura et al.l |2003| ). the global occurrence rate of su- 
pernovae accompanying magnetar remnants is -RsN.mag ~ 
IQ-'^Mpc-^yr-^ 

The intensity of the diffuse nonthermal neutrino back- 
ground is therefore roughly estimated as 

d$^ c dN^ 

—, ~ -. ; -rlSN,magtir?SF 

= 1 X 10-*C(e..o)CsFe;,o GeV-'cm-2 s"' sr"\(18) 

where tn = 10^" yr is the Hubble time, and we neglected the 
decrease of neutrino energy due to cosmic expansion (i.e., en- 
ergy redshift) for simplicity. We have optimistically assumed 
that each magnetar producing supernovae yields neutrino 
emissions comparable to our He star model. It is natural to 
assume that the magnetar birth rate, associated with the 
deaths of short-lived massive stars, is proportional to the 
cosmic star formation rate; the correction factor ^sf takes 
this into account. Since the star formation rate was larger 
by an orde r of magnitude in the pas t Universe at redshifts 
~l-2 (e.g., lHopkins fc Beacomll2006l), fsF might be about a 
few (see, for a similar discussion. I Waxman fc Bahcalllll999l ). 
We compare this intensity with that of atmo- 
spheric neutrinos. Aro und 0.1 GeV , the latter is ~ 
(ICaisser et all 1 19881 : iDaum et all 
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1 19951 : iMalek et all |2003| ). It is best to use the lowest- 
energy bin, around 0.1 GeV, to maximize the signal- 
to-noise ratio. At this energy, the intensity is ~3 x 
10~^5sF,o.5 GeV~^ cm~^ s~^ sr~^. This is many orders of 
magnitude smaller than the intensity of atmospheric neu- 
trinos, which makes it difficult to use the diffuse nonthermal 
neutrino background for extracting information on model 
parameters. We also note that this intensity at 0.1 GeV is 
smaller than the exponential tail of the intensi ty of diffuse 
thermal neutrino backgro und from supernovae (JAndo et al.l 
l2003l : lAndo fc SatdlJOoi l. 



4 DISCUSSION AND CONCLUSIONS 

In this paper we investigated proton acceleration and 
high-energy neutrino emission from the core-collapse su- 
pernovae of strongly magnetized stars. The stellar mag- 
netic field assumed in this work is determined under the 
fossil field hypothesis of magnetar magnetic fields. Since 
the C O core is generally op aque to high-energy neutrinos 
(e.g., iHoriuchi fc Anddl2008r ). we consider only high-energy 
neutrinos from the epoch when the supernova shock propa- 
gates through the stellar envelope. We considered three pro- 
genitor models and showed that protons may be acceler- 
ated up to ~10* GeV (CO star), ~10* GeV (He star), and 
~10^ GeV (BSG) respectively. These are above the pion pro- 
duction threshold, and we can expect high-energy neutrinos 
from pion decay. In all progenitors studied, the maximum 
proton energy is limited by proton-proton collisional cool- 
ing. 

The neutrino signal will be detectable above back- 
ground atmospheric neutrinos by the Super-Kamiokande de- 
tector, if a small fraction > 6 x 10~^ of the explosion energy 
is channelled into accelerated protons. Here we have quoted 
the results for a He star explosion (Type lb supernova) of 
Eex-p = 10^^ erg occurring at 10 kpc; this gives A^^ > 1 (see 
Eq. (HU). The required fractions are slightly larger for ex- 
plosions of CO stars and BSGs. These neutrinos, detected 
after thermal neutrinos but prior to the optical supernova, 
act as probes of large stellar magnetic fields. 

The neutrino signal can also be detected by km'^ de- 
tectors such as IceCube, which are sensitive to higher en- 
ergy (> 10 GeV) neutrinos. The predicted total event num- 
bers for IceCube deepcore are ~600 (CO star explosion), 
~6600 (He star explosion), and ~200 (BSG explosion) re- 
spectively. The strong suppression in a CO star explosion 
is due to strong pion cooling by inverse-Compton on elec- 
tron synchrotron photons. Thus the CO star result depends 
on the total energy of accelerated electrons, which is not 
well known. We have assumed a large fraction, noting that 
a smaller value would increase the total event number sig- 
nificantly. On the other hand, the strong suppression in a 
BSG star explosion is due to a low maximum proton energy. 

While the magnetar birth rate is conservatively esti- 
mated as 10~''yr~^ in the Galaxy, it could be larger by 
an order of magnitude because of its detection inefficiency. 
In addition, of t he ^10 magnetar candidates, a few lie 
within 3 kpc awav (jKothes et al.ll2002l : iHuUeman et al.ll2004l : 
ICaensler et al.l 120051 ). predicting a prolific neutrino signal. 
Both of these facts are positive aspects for testing this 
scenario by burst nonthermal neutrino detection. On the 



other hand, the diffuse nonthermal neutrino background 
from these sources is much smaller than conventional at- 
mospheric neutrinos. 

We have focused on the epoch of supernova shock prop- 
agation in the stellar envelope. As the shock continues 
propagating outwards, it eventually enters optically thin 
matter, i.e., it crosses the photosphere. At this point, ra- 
diation leaks forwards, and the shock makes a transition 
from radiation mediated to collisional (jEnsman fc Burrows! 
Il992[). The collisiona l phas e may be brief, as investigated 
by IWaxman fc Loea (|200lh . The latter authors find that 
the growth rate of electromagnetic instabilities within the 
shocked material is larger than the collisional rate, and the 
shock becomes coUisionless. E mission of TeV neutrino s have 
been predicted at this stage (jWaxman fc Loebll200ll ). Our 
scenario can be distinguished from this emission by lower 
neutrino energies, as well as the no-detection of simultane- 
ous 7-rays from neutral pion decay. 

We have for simplicity assumed a spherical supernova 
shock. Another simplification is our treatment of the mag- 
netic field. In reality, different shock and field geometries 
(parallel and perpendicular orientations) yield different ac- 
celeration efficiencies, and further treatments of both geome- 
tries and dynamics are required to address the full spectrum 
and luminosity curve of the neutrino signal. These are be- 
yond the scope of this paper, but we add that our simple 
assumption is in part due to the lack of realistic models of 
the magnetic field inside stars. 

Regardless of our approximations, we can show that any 
detection of the nonthermal neutrinos will occur after ther- 
mal neutrinos (if the case of SN 1987A is standard, then 
0(10) hours after), and last R-t/vs ~ seconds to hours, 
depending on the progenitor size. The neutrino signature 
from Type II supernovae, with larger progenitor radii, are 
expected to be longer in duration compared to those from 
smaller Type Ibc explosions. Detection will also be feasible 
for optically dark supernovae in the Galactic Centre or in 
molecular clouds, which are heavily obscured by dust. 
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